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ABSTRACT 

The determination of the brown dwarf binary fraction may contribute to the under¬ 
standing of the substellar formation mechanisms. Unresolved brown dwarf binaries 
may be revealed through their peculiar spectra or the discrepancy between optical 
and near-infrared spectral type classification. We obtained medium-resolution spectra 
of 22 brown dwarfs with these characteristics using the X-Shooter spectrograph at the 
VLT. We aimed to identify brown dwarf binary candidates, and to test if the BT-Settl 
2014 atmospheric models reproduce their observed spectra. To find binaries spanning 
the L-T boundary, we used spectral indices and compared the spectra of the selected 
candidates to single spectra and synthetic binary spectra. We used synthetic binary 
spectra with components of same spectral type to determine as well the sensitivity of 
the method to this class of binaries. We identified three candidates to be combination 
of L plus T brown dwarfs. We are not able to identify binaries with components of sim¬ 
ilar spectral type. In our sample, we measured minimum binary fraction of 9 . 1('13 q%. 
From the best fit of the BT-Settl models 2014 to the observed spectra, we derived 
the atmospheric parameters for the single objects. The BT-Settl models were able to 
reproduce the majority of the SEDs from our objects, and the variation of the equiv¬ 
alent width of the Rbl (794.8 nm) and Csl (852.0 nm) lines with the spectral type. 
Nonetheless, these models did not reproduce the evolution of the equivalent widths of 
the Nal (818.3 nm and 819.5 nm) and KI (1253 nm) lines with the spectral type. 
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1 INTRODUCTION 

Stars are believed to be born in large stellar nurseries, which 
they eventually leave to form the field population. A large 
number of stars remain in binary or hierarchical systems. 
Multiplicity provides constraints on fundamental parame¬ 
ters, such as dynamical masses, essential to test atmospheric 
and substellar formation models. It is well known that the 

* Based on observations of ESO, using VLT/ESO, under the pro¬ 
grams 084.C-1092(A), 085.C-0862(A) and 085.C-0862(B). 
f E-mail: manjavacas@iac.es; manjavacas@mpia.de 


binary fraction decreases when decreasing mass. This frac¬ 
tion decreases from 80%-60% for O and B stars, to 40% for 
the M dwarfs (Janson et al. 2012). The decreasing trend for 
binarity seems to extend to the substellar regime. For L- and 
T- brown dwarfs, the binary fraction is estimated at about 
20% (Gizis et al. 2003; Burgasser et al. 2003; Bouy et al. 
2003; Close et al. 2003; Burgasser et al. 2007; Luhman et al. 
2007; Goldman et al. 2008). 

Based on Monte Carlo simulations, Allen (2007) deter¬ 
mined that 98% of the brown dwarf binaries have separa¬ 
tions smaller than 20 AU. Burgasser et al. (2007) pointed 
out that the peak of the separation distribution of brown 
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dwarfs is ~ 3 AU, which is very close to the limit of the 
high resolution imaging surveys. Allen (2007) estimated that 
a fraction of ~6-7% of brown dwarf binary systems have not 
been detected yet, as a consequence of observational biases. 
For instance, Joergens (2008) searched for low-mass stars 
and brown dwarfs binaries in the Chamaeleon star forming 
region using the radial velocity method, and concluded that 
the percentage of brown dwarfs binary systems with separa¬ 
tions below 1 AU is less than ~10% in this star forming re¬ 
gion. Blake et al. (2010) monitored a sample of 59 ultra-cool 
dwarfs with radial velocity and determined that the binary 
frequency of low-mass unresolved systems is 2.5li g%. 

Spectroscopic data provide also important constraints 
to atmospheric models. These models allow us to disentan¬ 
gle the effect of varying effective temperature, gravity, and 
metallicity on the spectral features. Below a effective tem¬ 
perature of ~2600 K, models predict that clouds of iron and 
silicate grains begin to form, affecting the opacity (Lunine 
et al. 1986, Tsuji et al. 1996, Burrows & Sharp 1999, Fod¬ 
ders 1999, Marley 2000, Marley et al. 2002, Allard et al. 
2001). Self-consistent atmospheric models, such as the BT- 
Settl models (Allard et al. 2012a) and the Drift-PHOENIX 
models (Helling et al. 2008), use cloud models where the dust 
properties do not require the definition of any other addi¬ 
tional free parameters other than gravity, effective temper¬ 
ature and metallicity. Synthetic spectra for a specific set of 
atmospheric parameters can be compared to empirical spec¬ 
tra. For instance, these models have been tested on spectra 
of young late-type objects (late-type companions and free- 
floating objects; Bonnefoy et al. 2010; Witte et al. 2011; 
Patience et al. 2012; Bonnefoy et al. 2013; Manjavacas et al. 
2014). 

In this paper, we present X-Shooter optical and near 
infrared spectroscopy of 22 peculiar ultra-cool dwarfs, with 
spectral types between L3 and T7. We aim to hnd unre¬ 
solved brown dwarf binary systems. Our sample consists of 
objects with a different spectral classification in the optical 
and in the near-infrared or peculiar spectra in comparison 
with objects of the same spectral type. We want to con¬ 
tribute to the census of unresolved low mass dwarfs, and at 
the same time provide new constraints on the BT-Settl 2014 
atmospheric models. In Section 2 we describe the procedure 
for the selection of candidates in our sample, we explain 
how the observations were performed, and the data reduc¬ 
tion procedure. In Section 3 we hrst perform a search for 
L+T binaries in our sample, then perform a simulation to 
estimate the efficiency and false positive rate of the spectral 
htting method we apply. We estimate the sensitivity of our 
method to detect spectral binaries with the same subspec- 
tral type. We also compare our targets with trigonometric 
distances in a color magnitude diagram (CMD) with the L, 
L-T transition and T brown dwarfs published by Dupuy & 
Liu (2012). CMD allows us to discover unresolved binaries 
and young brown dwarfs. In Section 4 we discuss the proper¬ 
ties of the binary candidates selected in Section 3. In Section 
5 we update the binary sample of very low mass objects. In 
Section 6 we investigate how the BT-Settl 2014 atmospheric 
models reproduce our spectra over the optical and the near- 
infrared. We compare with the results from the literature, 
the equivalent width provided by the models. Finally, in 
Section 7 we summarize our results. 


2 THE SAMPLE, OBSERVATIONS AND DATA 
REDUCTION 

2.1 Sample selection 

We selected a sample of 22 brown dwarfs found in the liter¬ 
ature, with optical spectral types between L3 and T7, that 
have discrepant optical and near-infrared classihcation, or 
peculiar spectra. Optical subtypes are typically earlier than 
the near infrared subtypes. These objects are candidates to 
be unresolved binaries. We selected the brightest such ob¬ 
jects to ensure sufficient SNR in a reasonable integration 
time (J < 16). Furthermore, to calibrate our results and 
confirm the reliability of our method, we added some known 
brown dwarfs systems, LHS 102B (Golimowski et al. 2004), 
formed by a L4.5 plus a L4.5, and SDSS J042348.56-041403.4 
(Burgasser et al. 2005), formed by a L6±l and a T2±l. Our 
list of targets and their physical properties taken from the 
literature are compiled in Table 1. 


2.2 Observations and data reduction 

Our targets were observed using X-Shooter (Wideband 
ultraviolet-infrared single target spectrograph) on the Very 
Large Telescope (VLT) between October 2009 and June 
2010. X-Shooter is composed of three arms: UVB (300-550 
nm), optical (550-1000 nm) and near-infrared (1000-2500 
nm). It was operated in echelle slit nod mode, using the 1.6" 
slit width for the UVB arm, and the 1.5" slit width for the 
optical and the near-infrared arms. This setup provides res¬ 
olutions of ~3300 in the UVB and NIR, and ~5400 in the 
VIS. We obtained an average signal to noise of ~30. Ob¬ 
servations were performed at the parallactic angle to miti¬ 
gate the effect of differential chromatic refraction. We moved 
the object along the slit between two positions following an 
ABBA pattern with a size of 6 arcsec. The flux expected 
in the UVB arm is extremely low, therefore we chose not 
use spectra taken in this range. Telluric standards were ob¬ 
served before or after every target at a close airmass (±0.1 
with respect to the targets). Bias, darks and flats were taken 
every night. Arc frames were taken every second day. The 
observing log including telluric standard stars and the raw 
seeing during the observations is shown in Table 8. 

The spectra were reduced using the ESO X-Shooter 
pipeline version 1.3.7 (Vernet et al. 2011). In the reduction 
cascade, the pipeline deletes the non-linear pixels and sub¬ 
tracts bias in the optical or dark frames in the near-infrared. 
It generates a guess order from a format-check frame, a ref¬ 
erence list of arc line and a reference spectral format table. It 
rehnes the guess order table into an order table from an or¬ 
der definition frame obtained by illuminating the X-Shooter 
pinhole with a continuum lamp. The master flat frame and 
the order tables tracing the flat edges are created. Einally, 
the pipeline determines the instrumental response and sci¬ 
ence data are reduced in slit nodding mode. 

In the case of the near infrared, we used the spectrum of 
the telluric star of the corresponding science target observed 
in the same night to obtain the response function. We re¬ 
moved cosmetics and cosmic rays from the telluric stars, as 
well as the H and He absorption lines on their spectra, using 
a Legendre polynomial ht of the pseudo-continuum around 
the line. We then derived a response function by dividing 
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Table 1. List of observed targets: Magnitudes are in the 2MASS system, except for object G1 229B for which magnitudes are given by 
Leggett et al. (1999) in the UKIRT system. 


Number 

Name 

J [mag] 

H [mag] 

K [mag] 

dtTia (pc) 

SpT OPT 

SpT NIR 

Remarks 

Reference 

1 

LHS 102B 

13.11+0.02 

12.06+0.02 

11.39+0.02 

13.2+0.7 

L5 

L4.5 

Binary 

1, 2 

2 

2MASS J00361617+1821104 

12.47+0.02 

11.59+0.03 

11.31+0.02 

8.8+0.1 

L3.5 

L4 

NR“, 

3, 4, 36 

3 

2MASS J00531899-3631102 

14.45+0.02 

13.48+0.03 

12.94+0.02 


L3.5 

L4 


5, 6 

4 

SIMP 01365662+0933473 

13.46+0.03 

12.77+0.03 

12.56+0.02 

6.0+0.1 


T2.5 

V 

7, 8 

5 

2MASS J01443536-0716142 

14.19+0.02 

13.01+0.02 

12.27+0.02 


L5 


Red 

9, 10 

6 

2MASS J02182913-3133230 

14.73+0.04 

13.81+0.04 

13.15+0.04 


L3 

L5.5 


5,11 

7 

DENIS-P J0255.0-4700 

13.25+0.02 

12.20+0.02 

11.56+0.02 

4.9+0.1 

L8 

L9 

V 

12, 13, 35 

8 

2MASS J02572581-3105523 

14.67+0.03 

13.52+0.03 

12.88+0.03 

10.0+0.7 

L8 

L8.5 


4, 5, 14 

9 

2MASS J03480772-6022270 

15.32+0.05 

15.56+0.14 

15.60+0.02 

7.9+0.2 


T7 


15, 16 

10 

2MASS J03552337+1133437 

14.05+0.02 

12.53+0.03 

11.53+0.02 

9.1+0.1 

L5 

L3 

Y^ 

1, 17, 18, 19, 20, 37 

11 

SDSS J0423485-041403 

14.47+0.02 

13.46+0.03 

12.93+0.03 

15.2+0.4 

L7.5 

TO 

Binary 

1, 21, 33 

12 

2MASS J04390101-2353083 

14.41+0.02 

13.41+0.02 

12.82+0.02 

9.1+0.3 

L6.5 



11, 19 

13 

2MASS J04532647-1751543 

15.14+0.03 

14.06+0.03 

13.47+0.03 


L3pec 


Y? 

11, 14 

14 

2MASS J05002100+0330501 

13.67+0.02 

12.68+0.02 

12.06+0.02 


L4 

L4 


1, 22 

15 

2MASS J05395200-0059019 

14.03+0.03 

13.10+0.02 

12.53+0.02 

13.1+0.4 

L5 

L5 

NR 

1, 4, 24 

16 

2MASS J06244595-4521548 

14.48+0.02 

13.34+0.02 

12.59+0.02 

11.9+0.6 

L5pec 

L5 


1, 23 

17 

G1 229B 

13.97+0.03 

14.38+0.03 

14.55+0.03 

5.8+0.4 


T7pec 

MP*^, Y 

33, 34, 35 

18 

2MASS J10043929-3335189 

14.48+0.04 

13.49+0.04 

12.92+0.02 

17.0+1.6 

L4 

L5 


25, 26 

19 

2MASS J11263991-5003550 

14.00+0.03 

13.28+0.03 

12.83+0.03 


L4.5 

L6.5 

Blue L 

27, 28, 29 

20 

2MASS J13411160-3052505 

14.61+0.03 

13.72+0.03 

13.08+0.02 


L2pec 

L3 


22 

21 

2MASS J18283572-4849046 

15.18+0.05 

14.91+0.06 

15.18+0.14 

11.9+1.1 


T5.5 


23, 31 

22 

2MASS J21513839-4853542 

15.73+0.07 

15.17+0.09 

15.43+0.18 

16.7+1.1 


T4 


30 

References: [1] - Reid et al. (2008b), [2] - 

• Burgasser et al. (2007), [3] • 

- Dahn et al. 

(2002), [4] - Schneider et al. 

(2014), [5| - 

Marocco et al. 

(2013), [6] - 


Martin et al. (2010), [7] - Artigau et al. (2006), [8] - Radigan et al. (2013), [9] - Burgasser et al. (2011), [10] - Liebert et al. (2003), [11] - Liebert et al. 

(2003), [12] - Cruz et al. (2003), [13] - Castro et al. (2013), [14] - Kirkpatrick et al. (2008), [15] - Burgasser et al. (2003), [16] - Parker & Tinney (2013), [17] 

- Cruz et al. (2009), [18] - Allers Liu (2013), [19] - Faherty et al. (2013), [20] - Gagne et al. (2014), [21] - Antonova et al. (2013), [22] - Antonova et al. 

(2013), [23] - Faherty et al. (2012), [24] - Leggett et al. (2000), [25] - Andrei et al. (2011), [26] - Gizis (2002), [27] - Folkes et al. (2007), [28] - Faherty et al. 

(2009), [29] - Burgasser et al. (2008), [30] - Ellis et al. (2005), [31] - Burgasser et al. (2004), [32] - Vrba et al. (2004), [33] - Nakajima et al. (1995), [34] - 
Oppenheimer et al. (2001), [35] - Costa et al. (2006), [36] - Gelino et al. (2002), [37] - Zapatero Osorio et al. (2014). 

(a) NR: Not resolved binary; (b) V: Variability found; (c) Y: Young; (d) MP: Metal poor. 


the non-flux calibrated clean spectrum of the telluric stan¬ 
dard by a black body synthetic spectrum with the same 
temperature as the telluric star (Theodossiou & Danezis 
1991). Finally, to calibrate in response, we used the package 
noao.onedspec.telluric from the software Image Reduction 
and Analysis Facility, (IRAF). More details on data reduc¬ 
tion and flux calibration, as well as correction for telluric 
bands, are described in Alcala et al. (2014). 

To make sure that the flux in the whole near-infrared 
spectra was correctly scaled, we calibrated the flux of our 
near-infrared spectra using fluxes given by 2MASS (Two Mi¬ 
cron All Sky Survey). We convolved our near-infrared spec¬ 
tra with J, H and Ks Alter transmission curves of 2MASS. 
The resulting spectra were integrated. We calculated the 
flux for our targets corresponding to the J, H and Ks bands 
using 2MASS magnitudes (Cohen et al. 2003). Finally, we 
calculated the scaling factor for J, H and Ks bands and mul¬ 
tiplied our near-infrared spectra in J, H and Ks Alters to have 
the same flux as given by 2MASS. We scaled flux from the 
optical spectra to be consistent with the flux in the near- 
infrared. In the overlapping wavelengths of the optical and 
near-infrared spectra (995-1020 nm), we calculated a scaling 
factor, which is the median of the flux in the these wave¬ 
lengths of the near-infrared spectra, divided by the median 
of the flux in the overlapping wavelengths of the optical spec¬ 
tra. The reduced spectra are shown in Fig. 10^. Wavelengths 
affected by telluric absorption are removed from the figure, 
as well as the optical part for object G1229B, because it is 
contaminated by the flux of its companion, and the optical 
of 2M0144 because it is noisy. 


^ These spectra will be available in the ESO Phase 3 data release 


3 SEARCH FOR SPECTRAL BINARIES 

In this Section, we used different methods to reveal unre¬ 
solved brown dwarf binaries through their spectra. These 
methods are tailored to the type of brown dwarf binaries 
that we aim to And. 


3.1 Finding L plus T brown dwarf binaries 

The combined spectra of L plus T brown dwarf binary sys¬ 
tems are predicted to show peculiar characteristics. Those 
spectra are expected to have blended atomic and molecular 
absorptions of L and T brown dwarf spectra. This combi¬ 
nation may result in a peculiar spectrum. Burgasser et al. 
(2007, 2010) and Bardalez Gagliuffi et al. (2014) have stud¬ 
ied the spectral characteristics of L plus T brown dwarf bi¬ 
nary spectra and they have designed an empirical method 
to identify them using spectroscopy. 

Spectra of L plus T binary systems show bluer spectral 
energy distribution in the near-infrared than single objects 
of the same spectral type (Burgasser et al. 2010). Some spec¬ 
tral features vary: the CH 4 and H 2 O features at 1.1 /rm are 
deeper for binaries. The GH 4 feature at 1.6 /rm is stronger 
in comparison to the 2.2 /im CH 4 band. At 2.1 /rm the flux 
peak is shifted to the blue for the binaries. They also show 
larger flux from the T dwarf at 1.55 pm (Bardalez Gagliuffi 
et al. 2014). Using such differences, Burgasser et al. (2006, 
2010) and Bardalez Gagliuffi et al. (2014) defined spectral 
indices to identify L plus T brown dwarf binary candidates. 
Typically, spectral indices are defined as the ratio of spec¬ 
tral flux in two different wavelength intervals. The indices 
are specified in Table 9. Burgasser et al. (2006, 2010) and 
Bardalez Gagliuffi et al. (2014) compared all indices against 
each other for a large sample of brown dwarfs, some of them 
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L plus T known binaries. They identified the best pairs of in¬ 
dices that segregated known binaries from the rest of the ob¬ 
jects, and selected the regions in each combination of indices 
that delimit the known L plus T known binaries. In Table 10 
and Table 11 these regions are defined. There are several dif¬ 
ferences between Burgasser et al. (2006, 2010) and Bardalez 
Gagliuffi et al. (2014) methods. Burgasser et al. (2006, 2010) 
published eight indices that are valid only to find L plus T 
dwarf binaries, and defined six binary index selection crite¬ 
ria. The objects that satisfied two criteria were considered 
as "weak candidates". Those that satisfied three or more 
criteria are considered "strong candidates". Bardalez Gagli¬ 
uffi et al. (2014) used the eight spectral indices defined in 
Burgasser et al. (2006, 2010) and developed five new indices 
that are sensitive to M7-L7 plus T binaries. Bardalez Gagli¬ 
uffi et al. (2014) defined 12 new binary selection criteria. 
Objects satisfying four to eight criteria were considered as 
"weak candidates". Those that satisfied more than eight in¬ 
dices were considered "strong candidates". 

By calculating these spectral indices we selected those 
objects in our sample that are L plus T binary candidates. 
The result using Burgasser et al. (2006, 2010) criteria is 
shown in Figures 11, and the result using Bardalez Gagliuffi 
et al. (2014) criteria is shown in Figures 12 and 13. In Table 
2 , we summarize the weak and strong candidates given by 
each method. 

To confirm or reject the selected L plus T binary candi¬ 
dates, we compared our spectra with libraries of well charac¬ 
terized brown dwarf spectra, i.e. template spectra. We used 
as template spectra the McLean et al. (2003) and Cushing 
et al. (2005) libraries, with a resolution of R~ 2000, as well 
as the SpeX Prism Spectral Library spectra^, with a reso¬ 
lution of R ~120. In total we considered 462 spectra from 
SpeX Spectral Library plus 14 from Cushing et al. (2005) 
library and 47 spectra from McLean et al. (2003) library, 
with spectral types from LO to T7. 

We degraded the resolution of our X-Shooter spectra 
to the resolution of each template. We re-interpolated the 
library of brown dwarf template spectra and X-Shooter spec¬ 
tra to the same grid. We searched for the best matches to 
template spectra of single objects from SpeX, Cushing et al. 
(2005) and McLean et al. (2003) libraries and to synthetic 
binary spectra created using those libraries. To create those 
synthetic binaries, we calibrated the fluxes of the compo¬ 
nents to the same distance using an absolute magnitude- 
color relation (Dupuy & Liu 2012) and add them together. 
The final resolution of our synthetic binary templates was 
the same as the SpeX spectral library spectra. 

To identify the best matches to our spectra, we used 
the approach explained in Cushing et al. (2008), which is 
similar® to a 


® http://pono.ucsd.edu/~adam/browndwarfs/spexprism/ 

^ G is mathematically similar to a x^, but it does not follow a x® 
distribution as our comparison spectra have noise (see Cushing 
et al. (2008) for further details). We therefore do not expect to 
achieve G ~1 for the best fits, our goal is to determine whether 
a binary template is fitting better than a single template for our 
selected binary candidates. 



where C'(A) is the spectrum of the candidate, T(A) is the 
template spectrum, w(A) is a vector of weights proportional 
to the waveband size of each pixel, a is a scaling factor 
that minimizes G, and o'c(A) are the errors of the spectrum. 
To calculate the G, we used the parts of the spectra where 
no strong telluric absorptions are contributing, since we are 
confident of the telluric correction resulting from the data 
reduction process: A= 950-1350 nm, 1450-1800 nm and 2000- 
2350 nm. We additionally checked the best matches by visual 
inspection. Finally, we tested if the fit to a binary template 
was significantly better than the fit to a single template us¬ 
ing a a one-sided F-test statistic. We used as the distribution 
statistic ratio: 


single^ binary 

VSB = —-r-77- (2) 

TTllTlyijTifinary ) single 

where mm{Gsingie) and mm{Gbinary) are the minimum G 
for the best match to a single or to a composite template, 
and dfhinary and dfsingle are the degrees of freedom for the 
binary template fit and the single template fit. The degrees 
of freedom are the number of data points used in the fit 
(n = 296) minus one to account the scaling between our spec¬ 
tra and the template spectra. To rule out the null hypothe¬ 
sis, meaning that the candidate is not a binary with a 99% 
confidence level, we require rjsB > 1-31 The F-test anal¬ 
ysis rejected three of our candidates, namely: SIMP0136, 
2M0257, 2M1341. 

In Table 3 we show the best matches of the selected 
candidates to single and composite brown dwarf spectra. 
Plots with the best matches are shown in Appendix 7. 

We intended to estimate the fraction of missed L plus T 
binaries applying Burgasser et al. (2007, 2010) and Bardalez 
Gagliuffi et al. (2014) method. To this aim, we compared 
47 synthetic L plus T binaries to single L dwarfs, single T 
dwarfs and to other synthetic L plus T binaries. We found 
that 21% of the L plus T synthetic binaries, did not satisfy 
the binarity condition (i.e. they had r]sB < 1.31). In partic¬ 
ular, for 10 L plus T artificial binaries the match to single L 
dwarfs was significantly better than with binaries (see Fig. 
1). In Fig. 2, for all L plus T synthetic systems the bina¬ 
rity criteria was satisfied. Therefore, most of the L plus T 
systems should be found using this method, but it must be 
taken into account that some binaries might be lost. 

Finally, we compared a sample of 43 single L dwarfs to 
other L single dwarfs, and to synthetic L plus T binaries 
(see Fig. 3). We found that 37% of the single L dwarfs sat¬ 
isfied the binarity criteria, i.e. they had significantly better 
matches with L plus T synthetic binaries, they are therefore 
false positives. Equally, we performed a similar analysis for a 
sample of 40 single T dwarfs. We obtained that 35% of the T 
dwarfs are also false positives (see Fig. 4). We examined the 


Burgasser et al. (2006, 2010) requires a confidence level of 99%, 
and Bardalez Gagliuffi et al. (2014) a confidence level of 90%. We 
employed a confidence level of 99% to be more conservative, and 
minimize the false positives rate. 
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Table 2. Candidates selected by Burgasser et al. (2006, 2010) and Bardalez Gagliuffi et al. (2014) indices. 


Number 

Candidate 

Number of satisfied criteria from 
Burgasser et al. (2006, 2010) 

Type of candidate 

Number of satisfied criteria from 
Bardalez Gagliuffi et al. (2014) 

Type of candidate 

3 

2M0053 

2 

Weak candidate 

6 

Weak candidate 

4 

S1MP01365 

4 

Strong candidate 

8 

Strong candidate 

7 

DE0255 

2 

Weak candidate 

7 

Weak candidate 

8 

2M0257 

2 

Weak candidate 

6 

Weak candidate 

11 

SD0423 

2 

Weak candidate 

7 

Weak candidate 

20 

2M1341 

2 

Weak candidate 

8 

Strong candidate 


Table 3. Best matches to objects selected as binary candidates by spectral indices 


Candidate 

Single best match spectrum 

Composite best match spectrum 

VSB 

Fig- 

2M0053 

2MASS J17461199+5034036 (L5) 

Kelu-l (L3p) + SDSS J120602+281328 (T3) 

1.35 

14 

SIMP01365 

SDSS J152103.24+013142 (T2) 

DENIS-PJ225210-173013 (L7.5, bin) + SDSS J000013+255418 (T4.5) 

0.55 

15 

DE0255 

SDSS J085234.90+472035.0 (L9.5) 

SDSS J163030.53+434404.0 (L7) + SDSS J103931.35+325625.5 (Tl) 

3.42 

16 

2M0257 

SDSS J104409.43+042937.6 (L7) 

2MASS J0028208+224905 (L7) + SDSS J204749.61-071818.3 (TO) 

1.23 

17 

SD0423 

SDSS J105213.51+442255.7 (TO.5) 

2MASS J15150083+4847416 (L6) + SDSS J125453.90-012247.4 (T2) 

3.23 

18 

2M1341 

GJ1048B (LI) 

GJ1048B (LI) + 2MASS J1217110-031113 (T7.5) 

1.26 

19 



Figure 1. Histograms of r]gg for comparison of L plus T synthetic 
brown dwarf spectra to single L and other synthetic L+T brown 
dwarfs spectra. The dashed black line indicates r/gg = 1.31. 


spectral characteristics reported in the literature for the sub¬ 
samples of 16 L dwarfs and 14 T dwarfs with best matches 
to synthetic L plus T binaries. We found that 5 from the 16 
L dwarfs, and 2 of the 14 T dwarfs had either peculiar spec¬ 
tra or different spectral classification in the optical and the 
near infrared. These results suggest that this method is ef- 
hcient finding different spectral type binaries, but it should 
be applied with caution, as some false positives might be 
found. 



2 4 6 8 10 


Figure 2. Histograms of tfgg for comparison of L plus T synthetic 
brown dwarf spectra to single T spectra and other synthetic L+T 
spectra. The dashed black line indicates r;s 3 = 1.31. 


3.2 Finding equal spectral type brown dwarf 
binaries 

We aimed to find equal spectral type brown dwarf binaries 
comparing to spectral templates. To this purpose, we have 
chosen 101 LO to L9 presumed single brown dwarfs spectra 
from the SpeX library. Additionally, we have created 60 syn¬ 
thetic brown dwarf binaries in which both components have 
similar spectral type, i.e. the same spectral type, but dif¬ 
ferent spectral sub-types. To create those, we chose several 
presumed single L brown dwarfs from the SpeX sample, and 
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Table 4. Summary of the results obtained from Fig. 1, 2, 3, 4. 


Type of objects 

First comparison objects 

Second comparison objects 

Best matches 

Fig. 

Synthetic L+T binaries 

Single L dwarfs 

Synthetic L+T binaries 

21% to single L dwarfs (false negatives) 

79% to L+T synthetic binaries 

1 

Synthetic L+T binaries 

Single T dwarfs 

Synthetic L+T binaries 

0% to single T dwarfs 

100% to L+T synthetic binaries 

2 

Single L dwarfs 

Single L dwarfs 

Synthetic L+T binaries 

63% to single L dwarfs 

37% to L+T syntetic binaries (false positives) 

3 

Single T dwarfs 

Single T dwarfs 

Synthetic L+T binaries 

65% to single T dwarfs 

35% L+T synthetic binaries (false positives) 

4 



Figure 3. Histograms of T]gg for comparison of L single brown 
dwarf spectra to other L single dwarfs and synthetic L+T brown 
dwarfs spectra. The dashed black line indicates 7733 = 1 . 31 . 


Figure 4. Histograms of rjgg for comparison of T single brown 
dwarf spectra to other single T dwarfs and synthetic L+T brown 
dwarfs spectra. The dashed black line indicates 7733 = 1 . 31 . 


we created synthetic binary spectra combining single brown 
dwarf spectra following the same procedure as in Section 

3.1. 

We compared the 101 single L dwarfs and the 60 artifi¬ 
cial L plus L dwarf binaries to other L single SpeX spectra, 
and to other synthetic L binaries. We determined the best 
match of the 101 single L and the 60 synthetic L plus L 
binaries calculating the G parameter as in equation 1, and 
we equally decided the significance of the best match us¬ 
ing equation 2. We made a similar analysis for 56 T single 
dwarfs, and 74 T plus T synthetic binaries. 

We compared the best matches to single L and T dwarf 
spectra, and the best matches to synthetic binaries with 
similar spectral type (L plus L and T plus T, respectively). 
To this aim, we calculated the 7733 parameter for the 101 
single L dwarfs, 56 single T and for the 60 L plus L and 
63 T plus T artificial binaries created. In Figures 5 and 6 , 
we represent in a red histogram single brown dwarfs, and 
artificial dwarf binaries in a blue histogram with horizontal 
lines. 

We found that 49% of the 101 L single dwarfs, and 62% 


of the 60 synthetic L binaries, satisfied the binarity criteria. 
In Figure 5, we show that the distribution of 7733 for single 
L and for synthetic L plus L binaries is similar. 

Equally, we performed a similar simulation for T dwarfs. 
We found that 57% of the 56 T single dwarfs, and 86 % of the 
56 T plus T synthetic binaries satisfied the binarity criteria. 
In Figure 6 , we show the distributions of 7733 for single T 
and for synthetic T plus T binaries. 

For the cases mentioned before, the distribution of the 
7733 value is the same for single and synthetic binaries with 
the same spectral types, but different sub-spectral types. 
Therefore, it is impossible to distinguish between single L 
or T dwarfs and synthetic L or T dwarf binaries. Further¬ 
more, for both cases, the best fits are usually other synthetic 
binary spectra. Additional data, such as parallax measure¬ 
ments, high-resolution imaging or high resolution spectra 
are necessary in order to find these systems. 

3.3 Photometric search for brown dwarf binaries 

In our 22 object sample, distances for 15 objects are avail- 
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Table 5. Summary of the results obtained from Fig. 5 and 6. 


Type of objects 

First comparison objects 

Second comparison objects 

Best matches 

Fig. 

Single L dwarfs 

Single L dwarfs 

Synthetic L+L binaries 

51% single L dwarfs 

49% L+L syntetic binaries (false positives) 

5 

Synthetic L+L binaries 

Single L dwarfs 

Synthetic L+L binaries 

38% single L dwarfs (false negatives) 

62% L+L synthetic binaries 

5 

Single T dwarfs 

Single T dwarfs 

Synthetic T+T binaries 

43% single T dwarfs 

57% T+T synthetic binaries (false positives) 

6 

Synthetic T+T binaries 

Single T dwarfs 

Synthetic T+T binaries 

12% single T dwarfs (false negatives) 

86 % T+T synthetic binaries 

6 



Figure 5. Histograms of r]gg for comparisons of single L brown 
dwarfs (in red) and artificial L plus L binaries (blue with lines), 
to other single L dwarfs and other L plus L synthetic binaries. 
The dashed black line indicates r]gg=1.31. 

able in the literature, with a precision of 10% or better. In 
Figure 7 we present a CMD showing the J — K color in the 
MKO (Manna Kea Observatory) photometric system versus 
absolute magnitude in the J band. In this Figure, we plot 
all brown dwarfs with know parallaxes (Dupuy & Liu 2012), 
the 15 objects of our sample with known parallaxes, and 
the color-absolute magnitude relationship by Dupuy & Liu 
(2012). The two known binaries in our sample (LHS102B 
and SD0423, targets 1 and 11 respectively) stand out over 
objects with their same spectral types and other one, the 
young object 2M0355 (target 10) is much redder as objects 
of its same spectral type because of its youth. For the rest 
of the objects we cannot draw clear conclusions as there are 
no clear outliers. 


4 INDIVIDUAL CANDIDATES 

Six objects in our sample were selected as binary candidates 
by spectral indices. After fitting, three of our candidates 


Figure 6. Histograms of rigg for comparisons of single T brown 
dwarfs (in red) and artificial T plus T binaries (blue with lines), 
to other single T dwarfs and other T plus T synthetic binaries. 
The dashed black line indicates r]gg=1.31. 

were rejected due to the confidence level being lower than 
99%. This leads to a final number of three selected candi¬ 
dates. 

One of our selected candidates, the confirmed binary 
SD0423, was studied by Burgasser et al. (2005) and was 
used as a test of consistency of the spectroscopic method by 
Burgasser et al. (2010). We do not discuss it here. 

4.1 Rejected candidates 

We do not consider the following targets as a binary candi¬ 
date in the rest of the paper. 

4.1.1 SIMP 01365662+0933473 

SIMP 0136 was discovered by Artigau et al. (2006) and clas¬ 
sified as a T2.5. Goldman et al. (2008) searched for com¬ 
panions using NACO/VLT, reaching a sensitivity of 0.2" 
(1-40 AU), but no companions were found. Artigau et al. 
(2009) detected photometric variability in the J and K 
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J-K (mag) 


Figure 7. Color-magnitude diagram in the MKO system showing 
brown dwarfs with measured parallaxes from (Dupuy & Liu 2012), 
and its color-absolute magnitude relationship together with its 
plus minus one sigma curves (dotted line). Our targets are shown 
in black. Objects have the same numbers as in Table 1. 


the F-statistic. When we compared this target with objects 
of similar spectral type in the CMD, no overluminosity was 
found. This result agrees with the non-binarity scenario. 


4.1.3 2MASS J13411160-3052505 

The target 2M1341 was discovered by Reid et al. (2008a). 
Faherty et al. (2009) published a distance of 24±2 pc. Kirk¬ 
patrick et al. (2011) classihed it as a peculiar L2. Bardalez 
Gagliuffi et al. (2014) compared it to several SpeX templates 
and concluded that this object could be a L1.2T0.3 plus a 
T6.3±1.0. 

In Section 3.1, target 2M1341 was selected as a weak 
candidate L plus T binary by Burgasser et al. (2006, 2010) 
Indices, but it was selected as a strong candidate by Bardalez 
Gagliuffi et al. (2014) indices. This candidate was rejected 
by our conservative F-test criterion. However, a detailed in¬ 
spection of the observed spectroscopic matches (see Fig. 19). 
reveals no satisfactory reproduction of the near infrared fea¬ 
tures, thus leaving open the multiplicity of 2M1341. High 
resolution observations are required to disentangle its true 
nature. 


bands with a modulation of ~2.4 h and an amplitude of 
50 mmag. Radigan et al. (2012) calculated the amplitude of 
the variability for an object similar to SIMP 0136 (2MASS 
J21392676-I-0220226, T1.5). If this variability were produced 
by a companion, it would be much smaller than the variabil¬ 
ity obtained. Therefore, we do not expect that it is caused 
by a companion. Apai et al. (2013) explained it as a mixture 
of thick and thin patchy iron and silicate clouds covering the 
surface of the object. 

The object SIMP 0136 was selected as a brown dwarf bi¬ 
nary candidate, but it was rejected by a F-statistic analysis 
in Section 3.1. Spectral indices used in Section 3.1 are suit¬ 
able to select peculiar spectral characteristics that appear 
usually in binary L plus T brown dwarf spectra. However, 
if variability is produced by a partial coverage of thick and 
thin clouds in the brown dwarf atmosphere, similar peculiar 
spectral characteristics would appear in brown dwarf spec¬ 
tra. 

A preliminary parallax of 166.2 ± 2.9 mas for 
SIMP 0136 was obtained from the NPARSEC program (ESO 
program 186.C-0756, Smart et al. 2013). Using this paral¬ 
lax, we placed the object in a CMD together with other L, 
L-T transition and T brown dwarfs with parallaxes (Dupuy 
& Liu 2012), as shown in Figure 7. We compared SIMP 0136 
to objects of similar spectral type. We did not find signifi¬ 
cant overluminosity, expected in the case of late L and early 
T brown dwarf binaries. This result is compatible with the 
rejected binary hypothesis by the F-test. 

4.1.2 2MASS ,702572581-3105523 

The target 2M0257 was discovered by Reid et al. (2008b). 
Kirkpatrick et al. (2008) classihed it as a L8 in the optical. 
Marocco et al. (2013) measured its trigonometric parallax to 
be TT = 99.7T6.7 mas. It was selected by spectral indices as a 
weak brown dwarf binary candidate, and it was rejected by 


4.2 Selected candidates 
4.2.1 2MASS J00531899-3631102 

The object 2M0053 was discovered by Reid et al. (2008b). 
Kirkpatrick et al. (2008) classihed it as a L3.5 in the optical. 
It was selected as a weak brown dwarf binary candidate. We 
found a best match with a combination of a L3 dwarf and a 
T3 dwarf. There is not parallax measurement available for 
this target. 


4.2.2 DENIS-P .70255.0-4700 

The target DE0255 was discovered by Martin et al. (1999) 
and it was classihed as a peculiar L6. Koen et al. (2005) re¬ 
ported evidence of variability in different timescales (1.7 and 
5 h). Morales-Galderon et al. (2006) concluded that DE0255 
may vary with a 7.4 h period at 4.5 pm, but it does not 
at 8 pm. Costa et al. (2006) reported an absolute parallax 
of TT = 201.4± 3.9 mas. Burgasser et al. (2008) classihed it 
in the optical as a L8 and in the near infrared as a L9. Fi¬ 
nally, Reid et al. (2008a) searched for multiplicity for this 
target using high-resolution NICMOS NICl camera imag¬ 
ing on the Hubble Space Telescope, but found no evidence 
of multiplicity. 

In Section 3.1, DE0255 was selected as a L plus T weak 
binary candidate. We found a best match for DE0255 to a 
composite spectra of a L7 plus T1 spectra (see Fig. 16). 

Using the Costa et al. (2006) published parallax, we 
plot DE0255 in a CMD as above (see Fig. 7). In case of a 
late-L and early-T binary scenario, we expect to hnd about 
~0.5 mag overluminosity comparing with objects of similar 
spectral type on a CMD. In this case, no overluminosity was 
found, weakening the binarity hypothesis for object DE0255. 
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5 VERY LOW MASS BINARY FRACTION 

In our peculiar sample of 22 objects, we found three L+T 
binary candidates. One of them has been confirmed by other 
authors using high resolution imaging (SD0423, by Bur- 
gasser et al. 2005). Other two objects were selected as weak 
binary candidates (2M0053 and DE0255). The binarity hy¬ 
pothesis is weakened for object DE0255, due to the lack of 
overluminosity in the CMD (see Figure 7), expected for late 
L and early T brown dwarf binary scenario. 

This result allowed us to estimate the minimum and the 
maximum L-l-T binary fraction for our sample, and samples 
selected using our same criteria (see Section 2). The min¬ 
imum L-l-T binary fraction® is estimated at 4 . 5 ti' 4 % (^1^® 
only conhrmed L-l-T binary is SD0423 over the whole 22 
targets sample). The L-l-T maximum binary fraction for our 
sample is estimated at 18 . 614 ^ 3 "^% (the three L-l-T binary 
candidates over the whole 22 targets sample). These pairs 
would have a mass ratio of q ^ 0.5 for ages between 1 Gyr 
and 5 Gyr (expected for most of the objects in this study). 
Our work is not sensitive to smaller mass ratios. The de¬ 
rived range for the L+T binarity coincides with the frac¬ 
tion of late-M stars of the solar neighborhood that host T- 
type companions. As summarized by Burgasser et al. (2015), 
there are two late-M stars with T-type brown dwarf com¬ 
panions, among 14 dwarfs with spectral types between M7 
and M9.5, and at 10 pc from the Sun (~14%). Therefore, 
it appears that the L dwarf primaries have T-type compan¬ 
ions with a similar frequency to the late-M objects, despite 
the fact that the former primaries are expected to be less 
massive than the latter for typical field ages. 

Regarding binaries that include L+L and T+T pairs in 
our sample, there is just one confirmed L4.5+L4.5 system in 
our target list (LHS 102B). LHS 102B cannot be detected 
using the methods we employed in this paper, given the 
limitations of the spectroscopic technique. 

We thus determined the minimum fraction of L+L, 
T+T, and L+T pairs for our sample to be O.lta'o^o (the 
confirmed brown dwarf binaries, SD0423 and LHS102B over 
the whole sample of 22 objects). In spite of the peculiarity 
of our sample, this lower limit agrees with other values re¬ 
ported previously by different groups (Burgasser 2007; Gold¬ 
man et al. 2008; Bardalez Gagliufh et al. 2014; Burgasser 
et al. 2015). 


6 COMPARISON TO THE BT-SETTL 
ATMOSPHERIC MODELS 

The X-Shooter spectra presented in this paper provide the 
possibility to compare with the BT-Settl atmospheric mod¬ 
els in a wide range of wavelengths (550-2500 nm). We used 
13 objects of our total sample to test the BT-Settl mod¬ 
els 2014. We excluded brown dwarf binary candidates and 
spectra with low signal to noise to avoid false results. 

The BT-Settl models account for the formation and 
gravitational settling of dust grains for a effective temper¬ 
ature (Teff) below «2700 K in the photosphere of the ob- 

® The uncertainties of the binarity fraction for samples with less 
than 100 objects are calculated using the method explained in 
Burgasser et al. (2003). 


jects, following the approach described in Rossow (1978). 
The models include 180 types of condensates via their in¬ 
teraction with the gas phase chemistry, depleting the gas 
from their vapor phase counterparts. Fifty-five of these grain 
species are included in the radiative transfer calculations to 
the extent to which they have not settled from the cloud 
layer. Log-normal grain size distributions with a standard 
deviation of 1 are used, where the characteristic grain size 
for each layer is determined from the equilibrium size de¬ 
rived by the cloud model. The original timescales approach 
of the Rossow (1978) model has further been extended to 
account for nucleation as an additional timescale, which is 
defined by assuming a hxed seed formation rate motivated 
by studies of cosmic ray interactions with the Earth atmo¬ 
sphere (Tanaka 2005). 

The cloud model is implemented in the PHOENIX multi¬ 
purpose atmosphere code version 15.5 (Allard et al. 2001), 
which is used to compute the model atmospheres and to 
generate synthetic spectra. Convective energy transport and 
velocities are calculated using mixing length theory with a 
mixing length of 1 . 6 - 2.0 pressure scale heights, depending 
on surface gravity (log g), and overshoot is treated as an 
exponential velocity held with a scale height based on the 
RHD simulations of Ludwig et al. (2002, 2006) and Freytag 
et al. ( 2010 , 2012 ); an additional advective mixing term due 
to gravity waves is included as described in Freytag et al. 
(2010). All relevant molecular absorbers are treated with 
line-by-line opacities in direct opacity sampling as in Allard 
et al. (2003b); regarding this, the molecular line lists have 
been updated as follows: water-vapor (BT2, Barber et al. 
2006), vanadium oxide from Plez (2004, priv. comm.), TiO 
line list from Plez (2008) and collision-induced absorptions 
of H 2 (Abel et al. 2011). Non-equilibrium chemistry for CO, 
CH 4 , CO 2 , N 2 , and NH 3 is treated with height-dependent 
diffusivity also based on the RHD simulation results of Frey¬ 
tag et al. ( 2010 ). 

6.1 Comparison to synthetic spectra 

In this section, we compared X-Shooter optical and near- 
infrared spectra to predictions of the last version of the BT- 
Settl atmospheric models (Allard et al. 2003a, 2007, 2011) of 
2014. We exclude brown dwarf binary candidates and spec¬ 
tra with low signal to noise. We derive atmospheric param¬ 
eters of the objects and to reveal non-reproducibilities of 
the models. The models are described in Allard et al. (2011, 
2012 a,b). 

We selected subgrids of synthetic spectra with 400 
Teff SS 2100 K, 3.5 log g if 5.5 and metallicities of +0.0 and 
+0.3, which are the metallicities for which the latest version 
of the BT-Settl models are available. The solar metallicity is 
based on metallicities calculated by Caffau et al. (2011). The 
spacing of the model grid is 50 K and 0.5 dex in logg. Ef¬ 
fective temperature, gravity, metallicity and alpha element 
enhancement are described in the model name strings as 
lte-L0GG+[M/H]a+[ALPHA/H]. 

The BT-Settl 2014 synthetic spectra were smoothed to 
the resolution of X-Shooter. The models were then reinter¬ 
polated on the X-Shooter wavelength grid. The spectra were 
normalized using the same method as in Section 3 and ex¬ 
plained in Cushing et al. (2008). The results from the fit were 
always double checked visually. The atmospheric parameters 
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Table 6. Atmospheric parameters corresponding to the best 
fit spectra or synthetic fluxes for our targets. We give 
Teff/log s/[M/H]. 


Name 


Teff 

log 9 

[M/H] 

2MASS 

J00361617+1821104 

1800 

5.5 

+0.0 

2MASS 

J02182913-3133230 

1800 

5.0 

+0.3 

2MASS 

J03480772-6022270 

950 

5.0 

+0.0 

2MASS 

J03552337+1133437 

1700 

4.0 

+0.3 

2MASS 

J04390101-2353083 

1800 

5.0 

+0.3 

2MASS 

J04532647-1751543 

1750 

5.5 

+0.0 

2MASS 

J05002100+0330501 

1800 

5.0 

+0.3 

2MASS 

J05395200-0059019 

1800 

5.5 

+0.0 

2MASS 

J06244595-4521548 

1700 

4.5 

+0.3 

2MASS 

J10043929-3335189 

1800 

5.0 

+0.3 

2MASS 

J11263991-5003550 

1900 

5.5 

+0.0 

2MASS 

J18283572-4849046 

1100 

5.5 

+0.0 

2MASS 

J21513839-4853542 

1100 

5.0 

+0.0 


The grid size in Tefj is 50 K and in log g is 0.5 dex. Uncertain¬ 
ties in Teff and log g are the same as the grid sizes for both 
parameters. 


corresponding to the best fit models are reported in Table 
6 . The parameters Teff, log g, and [M/H] have uncertain¬ 
ties of 50 K, and 0.5 dex respectively. These errors corre¬ 
spond to the sampling of the atmospheric parameters of the 
model grids. We avoid the following objects to test models: 
binary candidates (2M1341, DE0255, SIMP0136, 2M0053, 
2M0257), known binaries (LHS102B, SD0423), noisy spec¬ 
tra (2M0144) or targets with known nearby objects that may 
contaminate the spectra, like in the case of G1229B. 

The CH 4 and the FeH molecules opacities are still in¬ 
complete in the new BT-Settl 2014 models. Methane line 
opacities are based on the semi-empirical list of Homeier 
et al. (2003), which is highly incomplete in the H band and 
only supplemented with a small set of room-temperature 
transitions for the Y and J bands. Iron hydride causes ab¬ 
sorption features through the E A — V* A system between 
650 and 1600 nm, but in addition to this Hargreaves et al. 
( 2010 ) identified significant opacity contributions from the 
A^n — A^n system, which is not yet included in the list 
of FeH lines available to PHOENIX. This explains that the H- 
band is not well reproduced for any of the L or T brown 
dwarf spectra, and also the J band in the case of T brown 
dwarfs. For three of the L brown dwarfs, the best match 
is found for log g= 5.5, with solar metallicity, four of the 
L brown dwarfs have best matches with log g= 5.0, but 
[M/H] = +0.3. 

Best matches to the BT-Setll models are shown in 
Fig. 20 and 21. There are two of the L brown dwarfs that 
have best matches with low gravity models, 2M0355 and 
2M0624. Object 2M0355 is known previously to be young 
(Cruz et al. 2009; Allers & Liu 2007; Zapatero Osorio et al. 
2014), so we expect gravity to be lower. The result given 
by the models is consistent with the literature. There are 
no references of youth for object 2M0624. Furthermore, the 
target does not show extremely redder J-K color in Figure 7 
or signihcantly weaker alkali lines on Figures 8 and 9, as low 
gravity objects do. Best matches to T type brown dwarfs are 
always solar metallicity models. The best match to object 
2M1828 is to a model with high gravity. 


6.2 Comparing predicted and observed equivalent 
widths 

We measured the equivalent width of a variety of alkali 
lines with sufficient signal-to-noise in our spectra and we 
compared those values to predictions of the BT-Settl 2014 
models. In the optical, we measured the equivalent width 
of the Rbl (794.8 nm), Nal (818.3 nm), Nal (819.5 nm) 
and Csl (852.0 nm). In the near infrared, we measured the 
KI (1253 nm) line. 

In Figure 8 , we plot the equivalent width of the alkali 
lines in the optical for the objects in our sample, for objects 
from Chiu et al. (2006), Colimowski et al. (2004), Knapp 
et al. (2004) and Lodieu et al. (2015) versus their spec¬ 
tral types. We previously degraded the spectral resolution 
of the spectra and the models to the X-Shooter resolution. 
In Figure 9, we plot equivalent widths of the KI line for 
our objects. We overplot field objects (Cushing et al. 2005; 
McLean et al. 2003), objects that belong to TW Hydrae 
Association (TWA), young companions (Allers & Liu 2007; 
Bonnefoy et al. 2014), young /?-dwarfs and 7 -dwarfs as a 
comparison (Allers & Liu 2007). We previously degraded 
the resolution of the observational spectra and the models 
to R~700, which is the lowest resolution of all the spectra 
for which we calculate KI (1253 nm) equivalent width. 

We overplot the equivalent width predicted by the BT- 
Settl models 2014 for those alkali lines. We transform pre¬ 
viously the effective temperature of the models to spectral 
types, using the empirical relation between spectral types 
and effective temperature relation published in Stephens 
et al. (2009). 

In the optical, the theory of cool atmospheres predicts 
the disappearance of the alkali elements in neutral form 
at temperatures below the L-T transition. This is a con- 
secuence of the depletion of the alkali atoms into molec¬ 
ular compounds, and the veiling by silicate clouds form¬ 
ing above the line-forming level, where the BT-Settl mod¬ 
els 2014 over estimate the dust scattering. The apparent 
strength of the Nal subordinate lines decreases with spec¬ 
tra types from the early Ls through the T dwarfs. The 
equivalent width of Csl (852.0 nm) increases from LO to 
L9 and it is maximum for the early T brown dwarfs, and 
it weakens progressively from the early to the late T brown 
dwarfs. The BT-Settl models reproduce the evolution of the 
equivalent width with the spectral type for Rbl (794.8 nm) 
and Csl (852.0 nm) lines, but underestimate the equivalent 
width of Nal (818.3 nm) and Nal (819.5 nm) lines, espe¬ 
cially in early to mid-L dwarfs. These elements do not par¬ 
ticipate directly in the sedimentation and dust formation, 
therefore the offset between the predicted and the observed 
equivalent widths could be due to uncertainties in the cloud 
model or the dust opacity contributing to pseudo-continuum 
that defines the equivalent width. 

In the near infrared (see Figure 9), the equivalent width 
of the KI (1253 nm) has two peaks at around L4 and T4, 
with a minimum at about L 8 . This might reflect that as 
for the potassium, we see different atmospheric layers for 
the various subtypes (Faherty et al. 2014). Object 2M0355 
(object 10 ) has weaker alkali lines in the optical and in the 
near infrared, as it is a young object (Faherty et al. 2012; 
Allers & Liu 2007; Zapatero Osorio et al. 2014). The BT- 
Settl models reproduce the weakening of the KI (1253 nm) 
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L2 L4 L6 L8 TO T2 T4 T6 T8 
Spectral type 


Figure 8. Equivalent widths of the detected alkali lines in the op¬ 
tical for our targets (red stars), binary candidates or known bina¬ 
ries from our sample (black stars) and for objects with equivalent 
width available in the literature (blue circles). The young object 
2M0355 (number 10) is marked as green star. These equivalent 
widths come from Chiu et al. (2006); Golimowski et al. (2004); 
Knapp et al. (2004) and Lodieu et al. (2015). We overplot with 
colored half circles joined by colored dashed lines the EW pre¬ 
dicted by the BT-Settl models 2014 for different gravities. The 
pink half filled circles and dashes lines correspond to log <; = 4, 
the blue ones correspond to log g = 4.5 and the grey ones to log g 
= 5.0. 


line for low gravity objects, but overestimates the equivalent 
width of this line for field objects. In the L-T transition, the 
KI (1253 nm) depletion or molecular blanketing is overesti¬ 
mated. 

In Table 7, we report the measured equivalent widths of 
the alkali lines for our sample. We do not report the equiv¬ 
alent widths for those alkali lines that were not detected in 
some of the targets. 



M4 M6 M8 LO L2 L4 L6 L8 TO T2 T4 T6 T8 
Spectral type 


Figure 9. Equivalent widths of the KI alkali line at 1253 nm of 
our objects (red stars), and binary candidates or known binaries 
of our sample (black stars), compared to the equivalent widths 
of field brown dwarf, young companions, young brown dwarfs 
(/3 and 7 dwarfs) and members of the TW Hydrae Association 
(TWA). The young object 2M0355 (number 10) is marked with a 
green star. We overplot with colored half circles joined by colored 
dashed lines the EW predicted by the BT-Settl models 2014 for 
differente gravities. We use the same colour code as for Figure 8 . 


7 CONCLUSIONS 

We observed and analyzed medium resolution VLT/X- 
Shooter spectra of 22 brown dwarfs with spectral types be¬ 
tween L3 and T7. Objects in our sample have peculiar spec¬ 
tral characteristics or different classifications in the optical 
and in the near infrared. Two of them are known binaries, 
that allow us to test our analysis. 

Using Burgasser et al. (2006, 2010) and Bardalez Gagli- 
ufli et al. (2014) empirical methods, we selected six objects 
as potential L plus T binary candidates: SIMP0136, SD0423, 
DE0255 and 2M1341, 2M0053 and 2M0257. We compared 
these six objects with single field brown dwarfs (McLean 
et al. 2003; Cushing et al. 2005 and Spex libraries) and syn¬ 
thetic binaries. We found the best matches using a statisti¬ 
cal analysis similar to the yf analysis. Objects SIMP0136, 
2M0257 and 2M1341 were discarded as candidates. The bi- 
narity hypothesis is weakened for object DE0255, due to 
their lack of over luminosity in the CMD, expected for late- 
L and early-T dwarf companions. 

We tested the efficiency of the method described in Sec¬ 
tion 3, and the possible proportion of false positives intro¬ 
duced aplying this method. To this aim, we compared single 
L and T dwarf spectra, and synthetic L plus T spectra, to 
other L and T single dwarf spectra, and to other synthetic 
L plus T binaries. We obtained that most of the L plus T 
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Table 7. Equivalent widths in nm for alkali lines measured in the optical and in the near infrared. 


Name 

Rbl (794.8 nm) 

Nal (818.3 nm) 

Nal (819.5 nm) 

Csl (852.0 nm) 

KI (1253 nm) 

LHS 102B 

1.19+0.01 

0.31+0.01 

0.38+0.01 

0.69+0.01 

0.79+0.01 

2MASS J00361617+1821104 

0.53+0.01 

0 .21+0.01 

0.24+0.01 

0.34+0.01 

0.93+0.01 

2MASS J00531899-3631102 

0.60+0.09 

0.16+0.02 

0 .22+0.02 

0.35+0.01 

0.85+0.01 

SIMP 01365662+0933473 
2MASS J01443536-0716142 

<0.3 

0.06+0.02 

<0.06 

0.65+0.02 

0.43+0.09 

0.83+0.01 

0.82+0.02 

2MASS J02182913-3133230 

<0.09 

0.17+0.02 

0.19+0.02 

0.31+0.01 

0.79+0.01 

DENIS-P J0255.0-4700 

<0.18 

0.13+0.01 

0.19+0.02 

1.04+0.02 

0.45+0.01 

2MASS J02572581-3105523 

<0.54 


0.15+0.05 

0.53+0.03 

0.54+0.01 

2MASS J03480772-6022270 

<0.46 

<0.28 

<0.013 

0.47+0.09 

0.42+0.01 

2MASS J03552337+1133437 

<0.03 

0.08+0.03 

0 .12+0.01 

0.29+0.01 

0.28+0.01 

SDSS J0423485-041403 

0.77+0.27 

0.09+0.02 

0.07+0.01 

0.59+0.02 

0.67+0.01 

2MASS J04390101-2353083 

<0.28 

0.17+0.03 

0.21+0.04 

0.57+0.02 

0.72+0.01 

2MASS J04532647-1751543 

<0.15 

0.21+0.05 

0.19+0.04 

0.25+0.03 

0.70+0.01 

2MASS J05002100+0330501 

0.64+0.31 

0 .21+0.01 

0.24+0.01 

0.37+0.01 

0.81+0.01 

2MASS J05395200-0059019 

0.83+0.03 

0.19+0.01 

0 .20+0.02 

0.54+0.01 

0.89+0.01 

2MASS J06244595-4521548 

G1 229B 

<0.61 

0.16+0.02 

0.19+0.03 

0.48+0.02 

0.70+0.01 

0.24+0.01 

2MASS J10043929-3335189 

0.75+0.30 

0.18+0.02 

0 .22+0.02 

0.40+0.02 

0.85+0.01 

2MASS J11263991-5003550 

0.75+0.21 

0 .22+0.01 

0.24+0.02 

0.45+0.02 

1.05+0.01 

2MASS J13411160-3052505 
2MASS J18283572-4849046 
2MASS J21513839-4853542 

0.37+0.28 

0.19+0.02 

0.25+0.01 

0 .21+0.01 

0.60+0.07 

0.55+0.16 

0.64+0.01 

0.71+0.01 

0.98+0.01 


synthetic binaries satisfy the binarity criteria. Nevertheless, 
21% of the L plus T synthetic binaries were not detected, i.e. 
21 % of those synthetic binaries are missed when we apply 
the method described in Section 3 using a 99% confidence 
level. 37% of the single L dwarfs, and 35% of the single 
T dwarfs satified the binarity criteria as well. The brown 
dwarf binary candidates found with this method should be 
confirmed using additional data. 

We examined the possibility of Ending equal spectral 
type brown dwarfs binaries. We compared single and syn¬ 
thetic binary spectra with the same subspectral type, to 
other single and synthetic binary spectra. For both cases, 
we obtained best matches with synthetic binary spectra for 
most of the cases. Therefore, we concluded that we are not 
able to find equal spectral type binary systems using this 
method. Additional data, such as parallax measurements, 
high-resolution imaging or high resolution spectra are nec¬ 
essary in order to find these systems. 

We re-calculated a lower limit for the very low mass 
binary fraction of 9.ll3 g% for our sample. We found that 
at least of L and T objects in our sample 

may be unresolved binaries with one L and one T possible 
members. This corresponds to a mass ratio of g Is 0.5 for an 
age of a few Gyr (expected for most investigated objects). 
This percentage agrees with previous results. 

BT-Settl models 2014 were able to reproduce the ma¬ 
jority of the SEDs of our objects in the optical and in the 
near infrared. Nonetheless, these models usually failed to 
reproduce the shape of the H-band, due to incomplete opac¬ 
ities for the FeH molecule in BT-Settl 2014 models. Best 
matches to models gave a range of effective temperatures 
between 950 K (T7) and 1900 K (L6.5), a range of gravities 
between 4.0 and 5.5. Some of the best matches corresponded 
to supersolar metallicity. 

We measured the equivalent width of alkali lines with 
good signal to noise (Nal, KI, Rbl and Csl) in the optical 
and in the near infrared spectra. We concluded that in the 


transition from L to T spectral types, the Nal doublet at 
818.3 nm and 819.5 nm in the optical is the first to disappear, 
while the other alkalines are present in the optical and near 
infrared in the whole L to T spectral types. We overploted 
the equivalent widths predicted by the BT-Settl models for 
those lines. The BT-Settl models reproduce the evolution of 
the equivalent width with the spectral type for the Rb I and 
Csl lines, and the weakening of the KI line for the early 
L with low gravity. Nevertheless, the models underestimate 
the equivalent width of the Nal, lines in the optical, and 
overestimate the equivalent width for the KI line for field 
objects. These elements do not participate directly in the 
sedimentation and dust formation. Therefore, the differences 
between models and observational equivalent widths may 
be due to uncertainties in the cloud model or in the dust 
opacities. 

The optical and near infrared spectra reported in this 
paper will serve as templates for future studies in any of 
these wavelengths. In the near future, the Gaia satellite will 
release high precision parallaxes of more than one billion 
of objects in the Milky Way, including hundreds of brown 
dwarfs. These parallaxes will allow us to detect the over¬ 
luminosity of brown dwarf binaries with respect to single 
brown dwarfs. 
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Table 8. Observing log: DIT is the integration time in each position of the slit, and NINT is the number of exposures. 


Name 

Date 

Arm 

DIT (s) 

NINT 

Seeing (”) 

Airmass 

Notes 

LHS102B 

October 16, 2009 

VIS/NIR 

290/300 

4/4 

1.0 

1.05 


Hip000349 

October 16, 2009 

VIS/NIR 

6/5 

1/1 

1.16 

1.31 

BOV Telluric Standard 

2M J0036+1821 

November 7, 2009 

VIS/NIR 

290/300 

4/4 

1.1 

1.4 


Hipll2022 

November 7, 2009 

VIS/NIR 

6/5 

1/1 

1.4 

1.5 

B2IV Telluric Standard 

2M J0053-3631 

October 16, 2009 

VIS/NIR 

290/300 

4/4 

1.44 

1.05 


Hip000349 

October 16, 2009 

VIS/NIR 

6/5 

1/1 

1.15 

1.01 

B9V Telluric Standard 

SIMP J0136+0933 

December 14, 2009 

VIS/NIR 

290/300 

4/4 

0.9 

1.2 


Hip021576 

December 14, 2009 

VIS/NIR 

6/5 

1/1 

1.1 

1.05 

B3V Telluric Standard 

2M J0144-0716 

December 14, 2009 

VIS/NIR 

290/300 

4/4 

1.05 

1.15 


Hip021576 

December 14, 2009 

VIS/NIR 

6/5 

1/1 

1.05 

1.15 

B6V Telluric Standard 

2M J0218-3133 

January 4, 2010 

VIS/NIR 

290/300 

4/4 

1.05 

1.15 


Hip009534 

January 4, 2010 

VIS/NIR 

6/5 

1/1 

1.1 

1.1 

B6V Telluric Standard 

DE J0255-4700 

October 17, 2009 

VIS/NIR 

290/300 

4/4 

2.2 

1.14 


Hip009549 

October 17, 2009 

VIS/NIR 

6/5 

1/1 

2.2 

1.2 

B6V Telluric Standard 

2M J0348-6022 

October 16, 2009 

VIS/NIR 

290/300 

12/12 

1.7 

1.3 


Hip012389 

October 16, 2009 

VIS/NIR 

6/5 

1/1 

1.7 

1.3 

B8V Telluric Standard 

2M J0355+1133 

December 21, 2009 

VIS/NIR 

290/300 

4/4 

0.92 

1.2 


Hip023060 

December 21, 2009 

VIS/NIR 

6/5 

1/1 

0.92 

1.2 

B2V Telluric Standard 

SD J0423-0414 

December 26, 2009 

VIS/NIR 

290/300 

4/4 

0.8 

1.4 


Hip020424 

December 26, 2009 

VIS/NIR 

6/5 

1/1 

0.9 

1.4 

BOV Telluric Standard 

2M J0439-2353 

December 21, 2009 

VIS/NIR 

290/300 

5/5 

1.4 

1.0 


Hip018926 

December 21, 2009 

VIS/NIR 

6/5 

1/1 

1.4 

1.0 

B3V Telluric Standard 

2M J0453-1751 

December 21, 2009 

VIS/NIR 

290/300 

8/8 

1.1 

1.1 


Hip023060 

December 21, 2009 

VIS/NIR 

6/5 

1/1 

1.1 

1.1 

B2V Telluric Standard 

2M J0500+0330 

February 05, 2010 

VIS/NIR 

290/300 

4/4 

0.7 

1.1 


Hip037623 

February 05, 2010 

VIS/NIR 

6/5 

1/1 

0.7 

1.1 

B5V Telluric Standard 

SD J0539-0059 

January 17, 2010 

VIS/NIR 

290/300 

4/4 

0.7 

1.1 


Hip033007 

January 17, 2010 

VIS/NIR 

6/5 

1/1 

0.7 

1.1 

B4V Telluric Standard 

G1229B 

December 14, 2009 

VIS/NIR 

290/300 

4/4 

1.2 

1.4 


Hip044786 

December 14, 2009 

VIS/NIR 

6/5 

1/1 

1.2 

1.4 

B6V Telluric Standard 

2M J0624-4521 

December 16, 2009 

VIS/NIR 

290/300 

5/5 

0.8 

1.4 


Hip030175 

December 16, 2009 

VIS/NIR 

6/5 

5/5 

0.8 

1.4 

B9.5V Telluric Standard 

2M J1004-3335 

February 5, 2010 

VIS/NIR 

290/300 

4/4 

0.9 

1.0 


Hip057861 

February 5, 2010 

VIS/NIR 

6/5 

1/1 

0.9 

1.0 

B5V Telluric Standard 

2M J1126-5003 

February 4, 2010 

VIS/NIR 

290/300 

4/4 

1.6 

1.1 


Hip073345 

February 4, 2010 

VIS/NIR 

6/5 

4/4 

1.6 

1.1 

B5V Telluric Standard 

2M J2151-4853 

May 08, 2010 

VIS/NIR 

290/300 

10/10 

1.4 

1.2 


Hiplll085 

May 08, 2010 

VIS/NIR 

6/5 

10/10 

1.4 

1.2 

BOV Telluric Standard 

2M J1341-3052 

June 02, 2010 

VIS/NIR 

290/300 

10/10 

0.7 

1.0 


Hip068124 

June 02, 2010 

VIS/NIR 

6/5 

10/10 

0.7 

1.0 

BOV Telluric Standard 

2M J1828-4849 

June 06, 2010 

VIS/NIR 

290/300 

10/10 

1.5 

1.4 


Hip092687 

June 06, 2010 

VIS/NIR 

6/5 

1/1 

1.5 

1.4 

B4III Telluric Standard 


APPENDIX 


OBSERVING LOG 
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REDUCED SPECTRA 

SPECTRAL INDICES CRITERIA 

BEST MATCHES TO POTENTIAL L PLUS T 
BINARIES 

BEST MATCHES TO BT-SETTL MODELS 2014 


This paper has been typeset from a LliX/ KrjX file prepared 
by the author. 
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Figure 10. Spectra of our 22 targets after reduction and degrading them at R~1000. Wavelengths largely affected by telluric absorption 
are removed from the figure in the near infrared, as well as the optical part for object G1229B, because it is contaminated by the flux of 
its companion and the optical part of 2M0144 because it is noisy. We plot spectra between 550-1350 nm, 1450-1800 nm and 1950-2500 nm 
to avoid telluric absortions. 
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Figure 11. Spectral index selection. Numbers 1-22 correspond to our objects. The boxes shown with dashed lines mark the areas where 
the selection criteria of Table 9 are valid. The red stars represent objects satisfying more than four such criteria. 
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Figure 12. Spectral index selection. 
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Figure 13. Spectral index selection. 
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Table 9. Spectral indices to select L plus T brown dwarf binary candidates. 


Index 

Numerator Range‘s 

Denominator Range“ 

Feature 

Reference 

HaO-J 

1140-1165 

1260-1285 

1150 nm H 2 O 

1 

CH4-J 

1315-1340 

1260-1285 

1320 nm CH 4 

1 

HaO-H 

1480-1520 

1560-1600 

1400 nm H 2 O 

1 

CH 4 -H 

1635-1675 

1560-1600 

1650 nm CH 4 

1 

HaO-K 

1975-1995 

2080-2100 

1900 nm H 2 O 

1 

CH 4 -K 

2215-2255 

2080-2120 

2200 nm CH 4 

1 

K/J 

2060-2100 

1250-1290 

J — K color 

1 

H-dip 

1610-1640 

1560-1590 + 1660-1690'’ 

1650 nm CH 4 

2 

K-slope 

2.06-2.10 

2.10-2.14 

K-band shape/ CIA H 2 

3 

J-slope 

1.27-1.30 

1.30-1.33 

1.28 ^m flux peak shape 

4 

J-curve 

1.04-1.07-^1.26-1.29'= 

1.14-1.17 

Curvature across J-band 

4 

H-bump 

1.54-1.57 

1.66-1.69 

Slope across H-band peak 

4 

H 2 O-Y 

1.04-1.07 

1.14-1.17 

1.15 p,mH 2 0 

4 

Derived NIR SpT 



Near infrared spectral type^ 

1 


(a) Wavelength range in nm over which flux density is integrated; (b) Denominator is the sum of the flux in the two wavelength 
ranges; (c) Numerator is the sum of the two ranges; (d) Near infrared spectral type derived using comparison to SpeX spectra. 
References: [1] - Burgasser et al. (2006), [2] - Burgasser et al. (2010), [3] - Burgasser et al. (2002), [4] - Bardalez Gagliuffi et al. 
(2014). 

Table 10. Index criteria for the selection of potential brown dwarf binary systems 


Abscissa 

Ordinate 

Inflection Points 

HaO-J 

HaO-K 

(0.325, 0 . 5 ),(0.65,0.7) 

CH 4 -H 

CH 4 -K 

(0.6,0.35),(1,0.775) 

CH 4 -H 

KjJ 

(0.65,0.25),(1,0.375) 

HaO-H 

H-dip 

(0.5,0.49),(0.875,0.49) 

Spex SpT 

HaO-J/HaO-H 

(L8.5,0.925),(T1.5,0.925),(T3,0.85) 

Spex SpT 

HaO-J/CH 4 -K 

(L8.5,0.625hT4.5,0.825) 


Table 11. Delimiters for selection regions of potencial brown dwarf binary systems 


Abscissa 

Ordinate 

Limits 


SpT 

CH 4 -H 

Best fit curve: y = -4.3xl0-'‘a;'^+0.0253x J- 0.7178 


HaO-J 

CH 4 -H 

Intersection of: -0.08x+1.09 and x = 0.90 


HaO-J 

H-bump 

Intersection of: y = 0.16x+0.806 and x = 0.90 


CH 4 -J 

CH 4 -H 

Intersection of: y = -0.56x + 1.41 and y = 1.04 


CH 4 -J 

H-bump 

Intersection of: y = l.OOx + 0.24, x = 0.74 and y = 

0.91 

CH 4 -H 

J-slope 

Intersection of: y = 1.250x -0.207, x = 1.03 and y = 

1.03 

CH 4 -H 

J-curve 

Best fit curve: y = 1.245x2 - 1.565x -f 2.224 


CH 4 -H 

H-bump 

Best fit curve: y = 1.36x2 - 4.26x + 3.89 


J-slope 

H-dip 

Intersection of y = 0.20x -f 0.27 and x = 1.03 


J-slope 

H-bump 

Intersection of: y = -2.75x + 3.84 and y = 0.91 


K-slope 

H 2 O-Y 

Best fit curve: y = 12.036x2 -20.000x +8.973 


J-curve 

H-bump 

Best fit curve: y = 0.269x2 - 1.326 + 2.479 
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Figure 14. Best matches for object 2M0053 to single (left plot) and composite spectra (right plot) . We show in black our X-Shooter 
spectra. In the plot on the left, the blue spectrum belongs to the best single match (2MASS 17461199+5034036, Reid et al. 2008a). In 
the plot on the right, we show in red the composite spectrum, in green the spectrum of the primary (Kelu-1, Ruiz et al. 1997) and in 
blue the spectra of the secondary (SDSS 120602.51+281328.7, from Chiu et al. 2006). rjgB = 1-35. The flux is F(A). 
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Figure 15. Best matches for object SIMP0136 (T2.5) to single (SDSS J152103.24+013142.7, from Knapp et al. 2004) and composite 
spectra (DENIS-PJ225210.73-173013, from Kendall et al. 2004 and SDSS J000013.54+255418.6, from Knapp et al. 2004). In black our 
smooth X-Shooter spectrum. Colors are the same as in 14. rjsB = 0.55. The flux is F(A). 
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Figure 16. Best matches for object DE0255 (L9) using single (SDSS J085234.90+472035.0, from Knapp et al. 2004) and composite 
spectra (SDSS J163030.53+434404.0, from Knapp et al. 2004 and SDSS J103931.35+325625.5, from Chiu et al. 2006). Colors are the 
same as in 14. r]SB = 3.42. The flux is F(A). 
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Figure 17. Best matches for object 2M0257 to single (SDSS J104409.43+042937.6, from Knapp et al. 2004) and composite spectra 
(2MASS J0028208+224905, from Cutri et al. 2003, and SDSS J204749.61-071818.3, from Knapp et al. 2004). Colors are the same as in 
14. = 1.23. The flux is F(A). 
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Figure 18. Best match for object SD0423 (TO) using single (SDSS J105213.51+442255.7, from Chiu et al. 2006) and composite spectra 
(2MASS J1515008+484742, from Wilson et al. 2003, and SDSS J125453.90-012247.4, from Leggett et al. 2000). Colors are the same as 
in 14. rjQB = 3.423.23. The flux is F(A). 
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Figure 19. Best matches for object 2M1341 (L2, peculiar) to single (GJ1048B, from Gizis et al. 2001) and composite spectra (GJ1048B, 
from Gizis et al. 2001, and 2MASS J1217110-031113, from Burgasser et al. 1999). Colors are the same as in 14. t)sb = 1-26. The flux is 
F(A). 
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Figure 20. Best matches to BT-Settl models 2014 found using equation 1, as in Section 3.1. Effective temperature, gravity, metallicity 
and alpha element enhancement are described in the model name strings as lte-L0GG+[M/H] a+[ALPHA/H] . The flux is F(A). 
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Figure 21. Best matches to BT-Settl models 2014 found using equation 1, as in Section 3.1. Effective temperature, gravity, metallicity 
and alpha element enhancement are described in the model name strings as lte-L0GG+[M/H] a+[ALPHA/H] . The flux is F(A). 

































